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Abstract 


The  development  and  ineorporation  of  the  latest  enhaneements  to  the  AVAST  eode 
are  deseribed.  The  purpose  of  this  work  was  to  make  the  modeling  of  the  physieal 
environment  more  realistic,  while  ensuring  that  the  code  runs  as  efficiently  as 
possible.  To  this  end  several  new  features  have  been  added.  These  include  modifying 
the  high-frequency  Kirchhoff  scattering  method  in  order  to  allow  for  at  least  one 
reflection  and  upgrading  the  existing  boundary  element  surface  panel  integration 
routines.  The  contract  also  addresses  the  need  to  investigate  the  high  frequency  target 
strength  of  Manta  shapes. 

Resume 


L’ elaboration  et  1’ integration  des  ameliorations  les  plus  recentes  apportees  au  logiciel  AVAST 
sont  decrites  dans  le  present  document.  Le  but  des  travaux  etait  de  modeliser  le  milieu 
physique  de  maniere  plus  realiste,  tout  en  veillant  a  ce  que  le  logiciel  soit  execute  le  plus 
efficacement  possible.  A  cette  fm,  de  nouveaux  attributs  y  ont  etc  ajoutes.  Ils  consistent  a 
modifier  la  methode  de  diffusion  de  Kirchhoff  pour  les  hautes  frequences  de  fafon  a  permettre 
au  mo  ins  une  reflexion  et  a  ameliorer  les  routines  actuelles  d’ integration  des  panneaux  de 
surface  par  elements  de  frontiere.  Le  contrat  porte  egalement  sur  I’etude  de  I’indice  de 
reflexion  a  haute  frequence  des  mines  Manta. 
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Executive  summary 


Introduction 

Work  by  DRDC  and  Martec  in  underwater/structural  acoustics  has  resulted  in  the 
development  of  a  series  of  eomputer  programs,  collectively  named  AVAST,  for  the 
numerieal  prediction  of  the  acoustic  radiation  and  scattering  from  floating  or 
submerged  elastie  structures  immersed  in  either  infinite,  half-spaee  or  finite-depth 
fluid  domains.  AVAST  eombines  both  the  finite  element  method  for  the  strueture  and 
the  boundary  integral  equation  teehnique  for  the  fluid.  The  finite  element  method  is 
used  to  predict  the  natural  frequeneies  and  related  mode  shapes  of  the  structure  in- 
vaeuo.  The  boundary  integral  equation  method  is  used  to  generate  a  system  of 
equations  relating  struetural  displaeements  to  fluid  aeoustie  pressures.  In  an  attempt  to 
make  the  modeling  of  sound  radiated  and  scattered  from  structures  more  realistic, 
several  enhancements  were  incorporated  into  the  existing  AVAST  suite.  These 
inelude  modifying  the  high-frequeney  Kirchhoff  scattering  method  to  allow  for  at 
least  one  refleetion  and  upgrading  the  existing  boundary  element  surfaee  panel 
integration  routines.  This  report  also  ineludes  a  summary  of  a  series  of  AVAST  elastic 
target  strength  analysis  eonducted  on  Manta-like  mine  models. 

Results 

The  AVAST  software  was  upgraded  to  allow  for  the  predietion  of  high  frequency 
target  strengths  of  underwater  targets  in  shallow  water,  to  inelude  the  effeets  of 
multiple  seatterers  and  seeondary  reflections,  and  to  allow  for  modelling  of  thin 
bodies.  A  realistic  examination  of  the  target  strength  of  a  Manta-like  mine  was  also 
performed  demonstrating  the  inereased  target  eeho  strength  at  frequencies  close  to  the 
struetural  resonant  frequencies  (when  submerged). 

Significance 

This  upgraded  version  of  AVAST  will  be  used  to  improve  the  aeeuraey  of  target  eeho 
strength  predietions  of  all  underwater  targets  and  the  results  from  the  mine 
investigation  will  be  used  to  eompare  with  other  measured  or  predicted  data  when 
available  to  assess  the  aeeuraey  of  the  AVAST  software  for  this  class  of  problems. 

Future  plans 

The  improved  AVAST  will  be  used  to  assess  the  aeoustie  target  strength  of  CF 
submarines  and  the  elastie  target  strength  of  other  mine-like  targets. 
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Sommaire 


Introduction 

Le  programme  de  recherche  en  collaboration  entre  RDDC  et  Martec  sur  I’acoustique  sous 
marine  (structures  immergees)  a  donne  lieu  a  I’elaboration  d’une  serie  de  logiciels,  regroupes 
sous  r  appellation  AVAST,  permettant  de  calculer  le  rayonnement  acoustique  et  la  diffusion 
du  rayonnement  par  des  structures  elastiques  flottantes  ou  immergees  dans  un  fluide  infmi,  ou 
fmie,  occupant  un  demi  espace  ou  un  volume  defmi.  Le  logiciel  AVAST  conjugue  la  methode 
des  elements  finis  (MET),  pour  la  structure,  et  la  methode  des  equations  integrales  de 
frontiere,  pour  le  fluide.  La  methode  des  elements  finis  (MEF)  sert  a  prevoir  les  frequences 
naturelles  et  les  formes  connexes  de  la  structure  in  vacuo.  La  methode  des  equations 
integrales  de  frontiere  sert  a  produire  un  systeme  d’ equations  qui  relie  les  deplacements 
structuraux  aux  pressions  acoustiques  du  fluide. 

Dans  le  but  de  rendre  la  modelisation  du  son  emis  et  diffuse  par  les  structures  plus  realiste,  on 
a  recemment  apporte  plusieurs  ameliorations  a  la  suite  AVAST.  Celles  ci  consistaient 
notamment  a  ameliorer  la  methode  de  diffusion  de  Kirchhoff  pour  les  hautes  frequences  afin 
de  tenir  compte  d’au  mo  ins  une  reflexion,  et  a  ameliorer  les  routines  actuelles  d’ integration 
des  panneaux  de  surface  par  elements  de  frontiere.  Le  rapport  comprend  egalement  le  resume 
d’une  serie  d’analyses  de  Tindice  de  reflexion  par  des  cibles  elastiques  realisees  sur  des 
modeles  semblables  aux  mines  Manta,  a  Taide  du  programme  AVAST. 

Resultats 

Le  logiciel  AVAST  a  etc  ameliore  afm  de  calculer  Tindice  de  reflexion  des  cibles  emettrices 
de  hautes  frequences  en  eaux  peu  profondes,  en  vue  de  tenir  compte  des  effets  des  diffusions 
multiples  et  des  reflexions  secondaires  et  de  permettre  la  modelisation  des  structures  minces. 
Un  examen  realiste  de  Tindice  de  reflexion  des  cibles  semblables  aux  mines  Manta  a 
egalement  etc  realise  dans  le  but  de  demontrer  que  Tamplitude  de  Techo  des  cibles  est  plus 
grande  a  des  frequences  proches  de  la  frequence  de  resonance  de  la  structure  (cibles 
immergees). 

Portee 

La  version  amelioree  du  programme  AVAST  permettra  d’ameliorer  Texactitude  des  calculs 
de  Tamplitude  de  Techo  de  toutes  les  cibles  sous  marines,  et  les  resultats  de  Tetude  sur  les 
mines  seront  utilises  pour  comparer  ces  previsions  avec  d’autres  donnees  mesurees  ou 
calculees  lorsqu’elles  seront  disponibles,  dans  le  but  d’evaluer  la  pertinence  de  Tutilisation  du 
logiciel  AVAST  pour  la  resolution  de  ce  type  de  problemes. 

Recherches  futures 

Le  programme  AVAST  ameliore  sera  utilise  pour  evaluer  Tindice  de  reflexion  acoustique  des 
sous  marins  des  FC  et  Tindice  de  reflexion  d’autres  cibles  elastiques  semblables  a  des  mines. 

Brennan,  D.P.  2007.  Further  Enhancements  to  the  High  Frequency  Target  Strength 
Prediction  Capabilities  of  AVAST.  DRDC  Atlantic  CR  2007-216.  DRDC  Atlantic. 
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1.  Introduction 


Phases  one  though  fourteen  of  the  DRDC/Martee  eollaborative  researeh  in 
underwater/  structural  acoustics  has  resulted  in  the  development  of  a  series  of 
computer  programs,  collectively  named  AVAST,  for  the  numerical  prediction  of  the 
acoustic  radiation  and  scattering  from  floating  or  submerged  elastic  structures 
immersed  in  either  infinite,  half-space  or  finite  depth  fluid  domains.  AVAST 
combines  both  the  finite  element  method  (FEM)  for  the  structure  and  the  boundary 
integral  equation  technique  for  the  fluid.  The  finite  element  method  is  used  to  predict 
the  natural  frequencies  and  related  mode  shapes  of  the  structure  in-vacuo.  The 
boundary  integral  equation  method  (BIEM)  is  used  to  generate  a  system  of  equations 
relating  structural  displacements  to  fluid  acoustic  pressures. 

In  an  attempt  to  make  the  modeling  of  sound  radiated  and  scattered  from  structures 
more  realistic,  several  enhancements  have  recently  been  incorporated  into  the 
previously  existing  AVAST  suite.  These  include  modifying  the  high-frequency 
Kirchhoff  scattering  method  in  order  to  allow  for  at  least  one  reflection  and  upgrading 
the  existing  boundary  element  surface  panel  integration  routines. 

In  the  discussion  which  follows,  details  concerning  the  development  and 
incorporation  of  these  latest  enhancements  to  the  AVAST  suite  will  be  presented. 

Also  included  in  this  report  is  a  summary  of  a  series  of  AVAST  elastic  target  strength 
analysis  conducted  on  Manta-shaped  mine  models. 
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2.  Manta  Shape  Target  Strength  Analysis 


The  basic  geometry  of  the  Manta  shape  is  displayed  in  Figure  2.1.  The  larger  radius  of 
the  shape  is  1000mm  and  the  smaller  radius  is  500mm.  The  shape  is  375mm  in  height. 
The  shape  was  modeled  using  4-node  quadratic  elements  and  the  mesh  used  is  shown 
in  Figure  2.2.  Figure  2.2  also  shows  the  interior  of  the  shape  with  the  core  section 
modeled.  The  outer  shell  of  the  shape  was  assumed  to  be  a  composite  material  with  a 
thickness  of  20mm.  The  core  and  the  bottom  plate  were  assumed  to  be  steel.  The  light 
blue  elements  represent  the  steel  portions  of  the  shape  and  the  brown  elements 
represent  the  composite  sections.  For  simplicity,  the  composite  was  modeled  as  an 
isotropic  material  with  a  density  of  2500kg/m^,  an  elastic  modulus  of  40GPa  and 
Poisson’s  ratio  of  0.35.  The  steel  elements  were  assigned  a  density  of  7900  kg/m^,  an 
elastic  modulus  of  207GPa  and  Poisson’s  ratio  of  0.30. 


2.1  Structural  Modeling 

The  algorithm  used  in  AVAST  for  computing  the  elastic  target  strength  is  based  on  a 
dry-mode  wet-mode  approach,  i.e.,  the  dry  (in-vacuo)  modes  of  the  structure  are  used 
to  predict  the  submerged  elastic  response.  For  the  purposes  of  this  analysis,  a  finite 
element  model  was  created  based  on  the  description  provided  above  and  a  series  of 
dry  structural  modes  were  computed  using  the  VAST  finite  element  code.  The  finite 
element  model  had  a  total  of  21084  degrees  of  freedom  with  a  total  of  30  modes 
computed.  The  results  of  the  natural  frequency  run  are  provided  below  in  Table  2.1. 
Note  that  for  the  purposes  of  this  analysis,  the  mine  was  modeled  as  a  free-free  system 
(no  boundary  constraints). 


2.2  Boundary  Element  Modeling 


The  boundary  element  model  used  to  represent  the  Manta  shape  is  represented  below 
in  Figure  2.1  (using  the  surface  panels  only)  and  represents  the  wetted  surface  of  the 
structure.  A  total  of  2552  quadrilateral  surface  panels  were  used.  Field  points  (for 
calculating  the  target  strength)  were  located  at  a  distance  of  100  meters  from  the 
center  of  the  shape  and  at  a  height  of  165  millimeters  above  the  sea  floor  (the  sea  floor 
was  represented  as  a  rigid  plane).  The  acoustic  source  used  to  excite  the  shape  was 
located  at  a  distance  of  1000  meters  and  a  height  of  165  mm  above  the  sea  floor.  The 
wet  natural  frequencies  computed  by  AVAST  are  provided  below  in  Table  2.2  and  the 
resulting  mode  shapes  are  shown  in  Figures  2.3-2.16. 
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Dry  Mine  Modes 

Wet  Mine  Modes 

Mode  Number 

Frequency  (Hz) 

Mode  Number 

Frequency  (Hz) 

7 

188 

7 

347 

8 

398 

8 

440 

10 

457 

10 

475 

11 

521 

13 

526 

13 

639 

15 

643 

15 

658 

17 

663 

17 

722 

19 

750 

19 

745 

25 

803 

21 

769 

27 

991 

23 

787 

29 

1005 

25 

795 

27 

904 

28 

983 

30 

1034 

2.3  Results 


In  order  to  study  the  effeets  of  struetural  elastieity  on  the  target  strength  of  the  Manta 
shape,  results  were  computed  for  narrow  8  Hz  frequency  bands  approximately 
centered  on  the  wet  natural  frequencies  of  the  structure.  Figures  2. 17  -  2.3 1  provide 
plots  of  the  radiated  component  of  the  acoustic  pressure  generated  by  a  unit  (1.0  Pa) 
incident  source  with  various  damping  factors  used  (from  0.0  to  0.5%).  The  horizontal 
axis  represents  frequency  while  the  vertical  axis  represents  field  point  position  in 
degrees  (one  complete  rotation).  These  plots  clearly  show  significant  changes  in  the 
radiated  pressure  at  the  frequency  approaches  a  coupled  resonance.  In  a  second  series 
of  figures  (see  Figures  2.32  -  2.44),  plots  representing  the  ratio  of  radiated  to  scattered 
acoustic  pressure  are  provided.  Once  again,  these  plots  clearly  show  a  significant 
increase  in  the  radiated  (i.e.,  elastic  contribution)  pressure  in  narrow  frequency  bands 
close  to  coupled  resonances. 


2.4  Target  Strength  Predictions 


In  order  to  compute  the  elastic  target  strength,  the  formula  used  to  compute  the  rigid 
target  strength  was  modified  to  include  the  radiated,  or  elastic,  acoustic  pressure,  i.e.. 


TS  Elastic  =  20Tog 


10 


^Rad  ^Scat 

<t>i 


+  lQLog,Q[R) 
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where  represents  the  seattering  field  pressure  represents  the  radiated  field 
pressure,  (j)-^  represents  the  ineident  pressure  eomputed  at  the  aeoustie  eenter  (i.e., 
eenter  of  the  mine),  and  R  represents  the  distanee  from  the  aeoustie  eenter  to  the  field 
points  (i.e.,  100  meters). 

Monostatie  elastie  target  strength  results  for  a  few  seleet  frequencies  are  provided 
below  in  Figures  2.45-2.50.  The  horizontal  axis  is  the  circumferential  direction  around 
the  mine  shape. 
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Figure  2.1 .  Schematic  Showing  Overaii  Dimensions  of  Manta  Shape 


Figure  2.2  .  Finite  Eiement  Mesh  of  Manta  Shape  with  Quarter  Cut  Out  to  Show  Core  Structure 
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Figure  2.3  .  Wet  Mode  7 


Mode  Shape  /Frequency  347.416 


Figure  2.4  .  Wet  Mode  7 
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Mode  Shape  8  Frequency  440.015 

Figure  2.5  .  Wet  Mode  8 


Mode  Shape  10  Frequency 


Figure  2.6 .  Wet  Mode  10 
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Mode  Shape  10  Frequency  474.774 


Figure  2.7 .  Wet  Mode  10 


Mode  Shape  13  Frequency  525.877 


Figure  2.8 .  Wet  Mode  13 
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Mode  Shape  15  Frequency 


Figure  2.9 .  Wet  Mode  15 


Mode  Shape  17  Frequency 


Figure  2.10 .  Wet  Mode  17 
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Mode  Shape  19  Frequency 


Figure  2.11 .  Wet  Mode  19 


Mode  Shape  25  Frequency 


Figure  2.12 .  Wet  Mode  25 
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Figure  2.13 .  Wet  Mode  25 


Mode  Shape  27  Frequency  991.224 


Figure  2.14 .  Wet  Mode  27 
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Figure  2.15 .  Wet  Mode  27 


Figure  2.16 .  Wet  Mode  29 
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Figure  2.17 .  Manta_365Hz_0e-3Damping 
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Figure  2. 18 .  Manta_379Hz_0e-3Damping 
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Figure  2.19 .  Manta_396Hz_0e-3Damping 
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Figure  2.20 .  Manta_419Hz_0e-3Damping 
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Figure  2.21 .  Manta_434Hz_0e-3Damping 
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Figure  2.22 .  Manta_434Hz_1e-3Damping 
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Figure  2.23 .  Manta_434Hz_2e-3Damping 
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Figure  2.25 .  Manta_532HzJ0e-3Damping 
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Figure  2.26 .  MantaJ545HzJ0e-3Damping 
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Figure  2.27 .  Manta Ja15Hz_0e-3Damping 
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Figure  2.28 .  Manta_660Hz_0e-3Damping 
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Figure  2.29 .  Manta_699Hz_0e-3Damping 
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Figure  2.30 .  iVianta_895Hz_0e-3Damping 
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Figure  2.31 .  Manta_900Hz_0e-3Damping 
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Figure  2.32 .  Manta_365Hz_0e-3Damping.RAT 
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Figure  2.33 .  Manta_379Hz_0e-3Damping.RAT 
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Figure  2.34 .  Manta_396Hz_0e-3Damping.RAT 


DRDC  Atlantic  CR  2007-216 


21 


Ratio  of  Radiated 
Over  Scattered 
Sound  Pressure 

3.858e+001 
3.390e+001 
2.921  e+001 
2.453e+001 
1.985e+001 
1.517e+001 
1.048e+001 
5.799e+000 
1.116e+000 


Figure  2.35 .  Manta_419Hz_0e-3Damping.RAT 
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Figure  2.36 .  Manta_434Hz_2e-3Damping.RAT 
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Figure  2.37 .  Manta_434Hz_5e-3Damping.RAT 
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Figure  2.38 .  Manta_532Hz_0e-3Damping.RAT 


DRDC  Atlantic  CR  2007-216 


23 


Ratio  of  Radiated 
Over  Scattered 
Sound  Pressure 


Figure  2.39 .  Manta_545Hz_0e-3Damping.RAT 
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Figure  2.40 .  Manta J615HzJ0e-3Damping.RAT 
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Figure  2.41 .  Manta_660Hz_0e-3Damping.RAT 
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Figure  2.42 .  Manta J699HzJ0e-3Damping.RAT 
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Figure  2.43 .  Manta_895Hz_0e-3Damping.RAT 
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Figure  2.44 .  Manta_900Hz_0e-3Damping.RAT 
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Target  Strength  (db)  Target  Strength  (db) 


Manta  Mine  Elastic  Target  Strength  at  379  Hz 


Figure  2.45 .  Manta  Shape  Monostatic  Eiastic  Target  Strength  at  379  Hz 


Manta  Mine  Eiastic  Target  Strength  at  419  Hz 


Figure  2.46 .  Manta  Shape  Monostatic  Eiastic  Target  Strength  at  419  Hz 
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Target  Strength  (dB)  Target  Strength  (db) 


Manta  Mine  Elastic  Target  Strength  at  615  Hz 


Figure  2.47 .  Manta  Shape  Monostatic  Eiastic  Target  Strength  at  616  Hz 


Manta  Mine  Eiastic  Target  Strength  at  660  Hz 


Figure  2.48 .  Manta  Shape  Monostatic  Eiastic  Target  Strength  at  660  Hz 
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Target  Strength  (dB)  Target  Strength  (db) 


Manta  Mine  Elastic  Target  Strength  at  699  Hz 


Figure  2.49 .  Manta  Shape  Monostatic  Eiastic  Target  Strength  at  699  Hz 


Manta  Mine  Elastic  Target  Strength  at  900  Hz 


Figure  2.50 .  Manta  Shape  Monostatic  Eiastic  Target  Strength  at  900  Hz 
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3.  Upgrade  the  AVAST  Boundary  Element  Surface 
Panel  Integration  Routines 


During  a  recent  series  of  numerical  trials,  it  was  found  that  the  numerical  integration 
algorithms  developed  for  use  in  AVAST  were  not  capable  of  maintaining  a  high  level 
of  accuracy  for  extremely  thin  bodies.  This  behavior  was  first  identified  in 
simulations  involving  ship  structures  having  highly  tapered  bows  (see  Figure  3.1  as  an 
example).  An  investigation  into  the  cause  of  this  loss  of  accuracy  determined  that  the 
numerical  integration  routines  used  to  generate  the  matrix  coefficients  used  in  the 
boundary  element  models  were  not  well-suited  for  evaluating  the  nearly-singular 
nature  of  these  coefficients.  This  nearly-singular  nature  occurs  when  the  centroids  of 
two  adjacent  panels  are  separated  by  a  relatively  small  distance  (relative  in 
comparison  to  the  dimensions  of  the  panels). 

In  order  to  correct  this  problem,  a  semi-analytic  approach  (based  on  the  Hess-Smith 
method  -  see  [1])  has  been  developed  and  implemented  into  the  latest  version  of  the 
AVAST  code.  This  new  approach  is  capable  of  computing  matrix  coefficients  for 
models  having  very  thin  or  even  collapsed  edges. 


Figure  3.1 .  Sample  Boundary  Element  Mesh  for  a  Ship  Hull 


30 


DRDC  Atlantic  CR  2007-216 


4.  Upgrade  Kirch h off  Scattering  Capabiiity 


During  the  previous  AVAST  development  eontraet,  a  version  of  the  Kirehhoff 
seattering  algorithm  was  implemented  in  the  AVAST  software.  While  target  strength 
predietions  generated  using  this  algorithm  appear  to  agree  elosely  with  other  eodes, 
the  eurrent  AVAST  formulation  does  not  aeeount  for  free  surfaee/shallow  water 
effects  or  reflections  from  adjacent  structure.  The  purpose  of  the  current  research 
work  is  to  address  these  modeling  deficiencies  and  produce  an  enhanced  version  of 
the  Kirehhoff  scattering  algorithm  in  an  upgraded  version  of  the  AVAST  code. 

Details  related  to  this  work  are  provided  in  the  discussion  which  follows. 

4.1  Free  Surface  -  Shallow  Water  Effects 

In  order  to  account  for  the  effects  of  free  surfaces,  rigid  surfaces  and  shallow  water 
fluid  domains,  a  formulation  based  on  the  method  of  images  [2]  has  been  incorporated 
into  the  latest  version  of  the  AVAST  code.  A  similar  formulation  has  been  employed 
in  AVAST  to  model  low  frequency  acoustic  radiation  and  scattering  via  the 
Helmholtz  equation.  In  cases  involving  free  or  rigid  surfaces,  a  single  “mirror”  image 
of  the  body  is  placed  on  the  opposite  side  of  the  boundary.  By  positioning  an  image 
in  this  way,  the  boundary  condition  on  the  boundary  (i.e.,  pressure  equals  zero,  or 
velocity  equals  zero)  is  automatically  satisfied.  A  similar  technique  is  used  for 
shallow  water  fluid  domains;  however,  a  series  of  images  must  be  used  (see  reference 

[3]). 

Fortunately,  only  minor  changes  were  needed  to  upgrade  the  AVAST  code  in  order  to 
provide  this  new  modeling  option.  In  fact,  the  bulk  of  the  effort  was  spent  modifying 
the  pre-existing  Kirehhoff  integral  formulation  (see  Equation  4.1)  in  order  to  account 
for  the  effects  of  the  secondary  images  (see  Equations  4. 2-4. 4). 


Infinite  fluid  domain: 

•  T  —ikRpQ 

Pscat  (^)  =  7^  f  - P:nc  scat  )  +  COs(^,„,  g  (4.1) 

4;r  J  RpQ 


Half  fluid  domain  with  free  surface: 

f  ~ikR 


p.Jp)  =  ^ls 


Q  -PQ  ^ 


v  ^PQ 


R 


Pine  {Ql^^A^scat  )  +  C0S(4^ 


PQ  ) 


(4.2) 
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Half  fluid  domain  with  rigid  surface: 
-ik 

■  +  ^ 


V  ^PQ 


R 


PinciOh^A^l^sca,)  +  COS 


(4.3) 


PQ  ) 


Shallow  water  (waveguide)  fluid  domain  of  depth  h\ 
-ik  ^ 

4^  -  V«=i 


_  jlr  (  ^  \ 

Pscat  (^)  =  ^  I  Z  P’  Pine  {Ql^^A^scat  )  +  COs(^,„,  )\lS q  (4.4) 

4^  i  Vtt 


where  P  represents  the  position  of  the  field  point,  Q  represents  the  location  of  the 
centroid  of  the  surface  panel,  k  represents  the  acoustic  wavenumber,  RpQ  represents 

the  distance  between  the  field  point  and  surface  panel,  RpQ  represents  the  distance 
between  the  field  point  and  image  surface  panel,  (q)  represents  the  magnitude  of 
the  incident  pressure  at  the  surface  panel  Q,  and  and  are  assumed  to 

represent  the  angle  between  the  surface  normal  and  the  vectors  representing  relative 
position  of  the  field  and  source  points  with  respect  to  panel  centroid. 


4.2  Multiple  Reflections 


In  order  to  account  for  multiple  reflections,  the  pre-existing  AVAST  Kirchhoff 
formulation  was  modified  so  that  the  scattering  pressure  computed  on  the  surface  of 
the  boundary  element  panels  would  be  modeled  as  secondary  source  terms.  In  terms 
of  the  AVAST  formulation,  given  the  particular  type  of  fluid  domain,  one  of  the  four 
equations  provided  above  would  be  used  to  compute  the  scattered  pressure 
contribution  made  by  source  at  the  various  field  locations.  This  contribution  is 
considered  the  “primary”  source  term.  Once  this  primary  contribution  has  been 
computed,  then  the  contributions  made  by  secondary  reflections  coming  from  the 
individual  panels  are  computed.  For  these  secondary  reflections,  the  location  of  the 
source  is  moved  to  centroids  of  the  individual  panels,  and  the  strength  of  the  source  is 
set  equal  to  the  scattered  pressure  computed  for  the  panel  centroid.  As  a  result,  for  an 
infinite  fluid  domain.  Equation  (4.1)  is  solved  once  using  the  strength  and  location  of 
the  primary  source  term,  and  then  repeated  once  for  each  surface  panel,  this  time 
using  the  location  of  the  panel  centroid  and  the  strength  of  the  scattered  pressured  at 
this  location  to  serve  as  the  secondary  source. 

4.3  Example 

In  order  to  demonstrate  the  capabilities  of  the  newly  upgraded  AVAST  Kirchhoff 
scattering  modeling  capability,  consider  the  model  of  boundary  element  model  of  a 
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generic  submarine  sail  provided  below  in  Figure  4.1  [4],  For  the  purposes  of  this 
study,  target  strength  predictions  have  been  made  for  two  model  configurations:  the 
sail  at  a  depth  of  3.0m  and  at  a  depth  of  300m.  For  the  300m  case,  the  fluid  domain 
was  assumed  to  be  infinite  in  extent.  Figure  4.2  provides  a  comparison  of  the  target 
strength  predictions  made  by  newly  upgraded  AVAST  code.  As  expected,  the  target 
strength  changes  significantly  (increasing  particularly  close  to  broadside  (90°))  as  the 
sail  approaches  the  free  surface. 


Figure  4. 1.  Generic  Submarine  Saii  Modei 
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Figure  4.2.  A  Comparison  of  Target  Strength  Predictions  for  Submarine  Saii  at  4  kHz  (infinite 
Fiuid  vs  a  Free  Surface  Height  3m  Above  the  Top  of  the  Saii) 
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